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NOMENCLATURE 

Cj anOlent speoU of sound* m/soc 

Cg group velocity* m/sec 

C} phase speed* m/sec 

f* frequency* Hz 

H unit step function, eqs. (22) and (23) 

1 nuflioer of axial grid points 

J numoer of transverse grid points 

Uesse) function of order tn 

L* length of duct, tn 

m spinning mode numoer 

P t1me>dependent dimensionless acoustic pres- 

sure, PiK,r,t), 

P^ time-dependent dimensionless acoustic pres- 

sure assoclateo with m mode 

Pm,n analytical solution for m spinning mode 
and n radial mode 

P* t Ime-depennent dimensionless acoustic pres- 

sure with angular venations, P'(x,r,«*t), 

P**/o5c52 
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ar 

T* 
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at 
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ax 

I* 

“m,n 
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«e 

c. 

’’cut 


nr 

0 
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spatially dependent "steady state* acoustic 
pressure, P(x,r) 

radial coordinate* r*/rQ 

radius of duct* m 

radial grid spacing 

period* l/f*, sec 

dimensionless time* t*/T* 

time step 

axial acoustic velocity, U*/Cq 
axial coordinate, x*/r5 
axial grid spacing 
acoustic Impedance* kg/m^ sec 
eigenvalue (n^h 2 «ro of Jm(a)) 
eq. (12) 
eq. (21) 

specific acoustic Impedance, 2*/o$/Co 

specific acoustic impedance at exit 

specific acoustic Impedance of Infinite duct 

cutoff freduency 

dimensionless frequency, r5f*/Cg 

angular coordinate 

dimensionless axial wavelength 


€ cutoff ratio, nr/ncut 

o5 aobtont air <ianatty* kg/mS 

w angular frequency 

Suoscripts: 

e exit condition 

1 axial Index (fig. 1) 

J radial index (fig. 1) 

m spinning node number 

n ('adial mode number 

0 ambient condition 
Superscripts: 

* dimensional quantity 

k time step 

INTROOUCVION 

Mith the introduction of strict aircraft noise 
regulations in tne late 1900* s the new aircraft 
nacelle designs required acoustic treatment in the 
inlet and exhaust ducts to reduce engine fan noise. 
To minimize the weight penalty of wall treatment, 
the aerospace industry ha*, been concerned with re- 
ducing the length of a liner for a required sound 
attenuation. Analytical techniques are needed to 
nandle sound propagation in ducts with axial varia- 
tions in cross-sectional area or in wall liner 
impedance (absorbers) and with gradients in the flow 
Mach number, in an attempt to meet this need, both 
finite difference and finite elemeni numerical tech- 
niques were oevelopeo. 

At the present time both "steady state" and 
transient numerical theories have oeen applied to 
the problem of sounu propagation in ducts. In the 
steady-state theory the pressure ana acoustic veloc- 
ities are assumed to be simple harmonic functions of 
time; thus tne equations governing sound propagation 
(linearized gas dynamic equations) becoiie indepeno- 
ent of time. Generally tne steaoy-state finite 
difference and finite element numerical algorithms 
nave been limited to low frequency and short ducts 
oocause of the large matrices associated with tne 
srlutions of the time-indepenoent equations. For a 
listing of recent publications and more details on 
tne techniques and proolems associateo with the 
numerical solutions of the sound propagation equa- 
tions, r»rer to reference 1, wnich contains a conv- 
prehensi e literature sumnary covering ooth finite 
difference and finite element analysis of small- 
ati^>litude (linear) sound propagation in straight a»u 
variable-area oucts. 

As an alternative to tne previously developed 
steady-state theories, time-depenaent numerical 
solutions were developed for plane-wave noise propa- 
gation III a two-dimensional rectangular duct without 

1 low (ref. ^), for parallel sheared mean flow 

(ref. i), and for higher order mode propagation in a 
cylindrical uuct witn a uniform mean flow (ref. 4). 
Advantageously matrix storage requirements are com. 
pletely eliminated in the time-dependent analysis. 
Unly tne solution vectors for pressure and velocity 


n«e<t be stored. Ttiys tlie transient solutions cen 
eeslly be applied In long ducts and In nifb- 
fraquency applications, tuclt as In tiw Inlet of e 
turbofan anginn* 

At the presmt t1«e, houever, the transient 
aethod appears to have one mjot draupack. Tne 
translant awthod, as fomulatod in reference 4, does 
not convergo for cutoff acoustic Mdos, This has 
Implications as to Its usa in a varlablt area where 
mooes may becoan cutoff In the small-erea portion of 
the duct, wnon a singla acoustic mode begins to 
move along a variable-area duct, part of its omirgy 
will be continuously transferred to the other duct 
modes (ref. b). If the initial mode is near cutoff, 
significant onargy could be transferred to a cutoff 
mode. Also, when abrupt changes in geometry or well 
impedance occur in a duct, a full ranga of acoustic 
mooes (including cutoff) is required to match pres- 
sure and velocity across the interface (ref. 6). 
Consequently the purpose of the present paper is to 
resolve the stability problem associated witn the 
transient calculation of cutoff acoustic modes. 

The numerical analysis to follow moduls the 
c losed-form pressure solutions for cutoff and propa- 
gating acoustic modes in a semi-ini .nite duct. Tne 
first section of the paper presents the various 
equations and boundary conditions governing sound 
propagation in a duct without a mean flow. Tne 
second sertion then presents tne difference form of 
the governiitg equations. Next analytical approxima- 
tions (steepest descent) are presented as a guide to 
tne understanding of the numerical results that 
follow. Numerical calculations are presented at 
forcing frequencies above, below, and nearly at tne 
cutoff frequency. At that time an explanation is 
presented for the computational instability asso- 
ciated witn cutoff modes. Finally ^ecoonendations 
are made on how to handle cutoff mooe propagation. 

GOVERNING EQUATIONS AND UOJNOARV CONDITIONS 

Tne propagation of sound in an axi iymroetric, 
cylindrical, naro-wall duct, as shown ii figure 1. 
is oescrioea by the wave equation and appropriate 
impedance boundary conditions. 

wave Equation 

The wave equation is a circular duct without a 
mean flow can be expressed in dimensionless form as 
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Tnese ano otner syinools are defined in tne nomencla- 
ture. The dimensionless frequency nr i& defined 
as 



Tne asterisks denote dimensional quantities. 

Oecause of the rotational nature or tne rotor 
blades on a typical turbofan Jet engine, large cir- 
cumferential variations in acoustic pressure will 
occur oepending on blade number and engine rpm. 

A three-dimensional solution for sound propagation, 
however, would be expensive to perform. Custom- 
arily, since tne equations are linear, tne cir- 
cumferential acoustic pressure variations are 
decomposed into spinning mooes m: 
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The sunudi ton U over those nodes thet ere present 
tn a particular application. Considering solutions 
with a single spinning mode number m 

P’U.r.e.t) • P(K,r,t)e*"* (4) 

Che wave equation (1) reduces to 
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Equation (b) in difference form will be solved to 
determine the pressure in the duct 

Haro-wail Boundary Condition 

The boundary condition at the surface of a 
hard-wall duct is 


ar 
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(6) 


Entrance Condition 

tne boundary condition at tne source plane 
P(0,r,t) will be assumed to vary as or in 

a dimensionless form as, furthermore the 

transverse pressure variation (r direction) will be 
assumed to correspond to the eigenfunctions 
dm(<*mn'') associated with mode propagation in an 
infinitely long nard-wall duct. The eigenvalues 
associated witn mode m,n are tabulated in refer- 
ence ? (p. bll) and reference 8 (p. 411). Therefore 
the source boundary condition used herein is 

PlO.r.t) . V<.^r)e‘^"‘ (7) 


Exit Impedance 

The bounoary condition at the exit of the duct 
can be expressed in terms of a specific acoustic 
impedance defined as 


«e • U 

where the equation describing tne acoustic velocity 
0 is 
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Substituting equation (8) into equation (9) yields 
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Ih reference 4, has been assumed to be 
the sieaoy-state impedance associated with mode 
propagation oown an infinite duct, for transmission 
of a single acoustic mode witnout reflection, tne 
exit impedance is 
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for propagating modes and 


for cutoff acoustic modes. The negative root is 

required for damping, when e^^***^* is assumed to 
be the forcing runction. 

All the other equations and boundary conditions 
are exact. The exit impedance was the only 
condition that was modeled in an approximate man- 
ner. As indicated by its title, this paper con- 
centrates on exact treatment of tne exit impedance 
condition. 

Centerline Condition 

for plane-wave propagation (m ■ 0, n »> 0), 
aP/ar is zero at t»d centerline; for modes with m 
equal to or greater than 1, P(x,0,t) is zero. Also 
aP/ar is zero for m greater than 1. 

Initial Condition 

for times equal to or less than zero the duct 
is assumed to be quiescent; that is, the acoustic 
pressures and velocities ar>7 taken to be zero, for 
times greater than zero tne application of tne noise 
source (eq. (7)) will drive the pressure in the duct. 

UlffEAENCE Ei)UATlONS 

Instead of a continuous solution in space and 
time tne finite-difference approximations will 
uetermine tne pressure at isolated grid points in 
space as shown in figure 1 and at discrete time 
steps at. Starting from the known initial condi- 
tions at t ■ 0 ana the boundary conditions, the 
finite-difference algorithm will march out tne solu- 
tion to later times. No special starting equation 
will be required, because both tne pressure and 
acoustic velocities are initially ass*imed to be 
zero, a quiescent duct. 

Oifference Equations 

Away fr^ the duct boundaries, in cell 1 of 
figure 1. the first and second derivatives in the 
wave equation (eq. (a)) can be represented by the 
usual central differences in time ana space (ref. 9, 
p. 4b2) 
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(14) can be rewritten as 
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Equation (IS) is an algorithm tnat permits marching 
out solutions from tcnoitn vahes of pressure at times 
associated wun k and k - 1. ilfference 
equations for cells ^ S in f > -, )r> 1 are cf 
sligntly different form. Tnese ,.at1ons as «call as 
tneir derivations are given In references I and 4. 

Spatial rtesh Size 

Vhe mesn spacings ax and ar must be re- 
stricted to small values to reduce the truncation 
error. To resolve the oscillatory nature of the 
acoustic pressure, tne required number of grid 
points I In the axial direction and J In tne r 
direction were given in reference 4 as 
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In tne explicit time-marching approach used 
nere roundoff errors can grow in an unoounoed fash- 
ion and destroy the solution If the tiine increment 
at IS too large or if tne iteration scneme is im- 
properly posed. For spinning wave propagation in a 
tnin annulus tne Von Neumann metnoa (ref. lO) 
applied to equation (14) yields 

0 - ar 

at < — ' " T' r y (Id) 
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For circular oucts at was empirically decreased oy 
a factor of Q.6. 

AIWM.VT1CAL SOLUTION 

Tne wave equation (1) and its associated bound- 
ary ana initial conditions can oe solved directly 
for tne semi-infinite duct shown in figure ^(al by 
the use of tne Laplace transform (ref. »). To oi> 
tain a more useful solution, Pearson (ref. 6) has 
also solved equation (1) oy tne metnod of steepe'.t 
decent to obtain 


P(M.r.t) 4.1.. _ _4 ftrtm1ent\ 
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where H Is the unit step function and 



The functional form of the transient solution Is 
presented later. 

In reference b tne duct eigenfunctions were 
left in a general form. Herein the cylindrical form 
of tne eigenfunction Is used, and the original 
equations of reference b are rewritten In a dimen- 
sionless form. Also, only the real part of the 
transient and steady-state solutions was presents^ 
in reference S. These solutions, however, are pre- 
sented in complex form in equations (14) and (20). 
decause tne transient terms correspond to waves that 
are propagating but are not simple narmonic in 
nature, the usual restriction that rufoff waves do 
not carry energy (ref. 11, p. 5J) not apply 
during the transient. 

As a result of the method of steepest decent 
the acoustic pressure field inside tne duct Is split 
into a transient motion and a steady-state solution 
that IS either propagated (eq. (14)) or attenuated 
(eq. (20)} depending on wnether the frequency of the 
forcing field is greater or lets than tne cutoff 
frequency of the mode considered. The steady-state 
solution is also the soUtion to equation (1) when 
tne pressure is assumed ‘.o be a simole narmonic 
function of time. 

For propagating modis the physical significance 
of the unit step functions in equation (14) can best 
be described if they are rewritten in dimensional 
form: 
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Tiuit tiM transimt tam In aquation (19) travals at 
taa tlqnal valoefty (ftf. 7. p. 479) utiKh It tha 
tpaa4 of found of ttia Mdlua. Tlw trmtlant front 
IS dafiMd at tlM turfact bayond Miica tha Mdlwi It 
coaplataly at ratt at any pivon Instant of tiai. On 
tha othar hand tha staa^stata solution travalt at 
tha group valocity Cg« uhlch Is simply 
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tiaos k and k - 1 Mod to bo storad tinea only 
tha prottwrat at M tmo pravlout tiaM ttaps art 
naadad for tha p^* calculation (aq. (IS)). 

Aftar sufflelant tloa has patsad to that tha staadiy* 
stats solution hat travalad to tha position x ■ 1. 
(H(t - nfjt) or H(t - nrly) ■ 1) and aftar an 
additional tint parlod aiaptat such that tha trait* 
slant taros hava diad out (chackad nuoarlcally). the 
tlon-depandent results can be compared with tha 
staadystata results. In this comparison tN time* 
dapandant pressure Is simply divided by a'<*' to 
obtain 


ttsf. 7^..'}.sr<iShS3 »<.. 

tions. Tha phase velocity C} Is always 
reator than Cq; the group velocity Is always 
ess than Cq. 

SIMULATION Of SCNI-INFINITE OUCT 

rne form of the transient solution leads to 
some observations that are pertinent to the numeri- 
cal results to follow. First, tha time In which to 
perform the transient solution Increases near the 
cutoff frequency. In fact, equation (23) Indicates 
that a mooe exactly at cutoff can never exist at 
steady state since its group velocity 1$ tero. 
Second, the Impedance at any position will not be 
equal to the steady-state Instance during the Ini- 
tial transient. Thus equation (11) Is not the exact 
analytical condition for no reflections during the 
transient. As a consequence of using equation (11) 
during the transient, part of the acoustic energy 
reaching the exit will be reflected back toward the 
source. 

For modes with small group velocities near the 
cutoff frequency, the possibility exists that the 
use of equation (11) as a boundary condition could 
Introduce an unstaole feedback loop Into tne numeri- 
cal analysis. This computational Instability could 
prevent tne estaoilshment of the true steady-state 
solution. Therefore tne possibility that the exit 
impedance was the cause of tne instability of cutoff 
acoustic modes (ref. 4), as discussed In the Intro- 
duction, is now Investigated. 

The present analysis attempts to numerically 
reproduce tne steady-state analytical solutions 
given by equations (IP) and (20) for the acoustic 
pressure in tne oomai'i 0 < x < 1 in tne semi- 
infinite duct of figu' e 2Ca).“ A variety of fre- 
quencies above and below the cutoff frequency are 
chosen. In tne numerical slnxilatlon a finite-length 
uuct (L*/rQ * 1), shown In figure 2(b), and a 
long ouct (L*/rg > 40), »;iown In figure 2(c), 
are used. The latter case exactly slmulates-an In- 
finite duct if the numerical calculations are ter- 
minated before the return of the Initial reflected 
wave from the duct exit. 

OLSCUSSION OF RESULTS 

In this section the calculated axial pressure 
profile of the 3,0 mode Is examined at frequencies 
above, below, and nearly at the cutoff frequency, 
in these problems the time-dependent numerical 
results are compared with the exact solutions given 
by tne second term In equations (IP) and (20). The 
staolllty of tne numerical solutions Is examined for 
source frequencies near tne cutoff frequency. 

In tne numerical computer program the values of 
pressure are calculated at each grid point l.j for 
each time step k. Only the values of pressure at 


p . (25) 

ft 

which represents a mimericel approximation to the 
Fourier transform of P(x,r,t) as stated In refer- 
ence 12 (p. 11). 

in all tne cases to be examined now the numeri- 
cal and analytical pressure profiles are displayed 
for the domain 0 < x < 1 for a hard-wall, semi- 
infinite duct witi) an~ ow3 spinning mode and the 
lowest order radial eigenvalue aj a * 4.20119). 

The values of the acoustic pressim are displayed 
at r ■ 1, along the duct wall. Also, to keep the 
graphical displays sliaple, only the real components 
of pressure are compart. 

Tne ratio of the forcing frequency nr To 
tne cutoff frequency ncut defined as e. 

For he 0^3 spinning mode, the cutoff frequency Is 
ncut " 0.66864. Therefore 

« ■ ofesT 

Finally, In the discussion to follow, the 
ratio of the speed of sound to the group velocity 

S /Cfl Is cited. This ratio Is equal to the 
melislonless time t required for the steady- 
state solution to propagate to the x-1 position in 
the duct. This parameter Is also used as a measure 
of how close the forcing frequency Is to the cutoff 
frequency. The difference In the calculated time 
t and the ratio C$/C« represents the added 
computational time required for the transient terms 
to die out so that the steady-state terms become 
dominant. 

Propagating Acoustic Modes 

The numerical and analytical profiles are com- 
pared in figure 3 for nr > 1. In this case the 
ratio of sonic to group velocity Cg/Cg 
is 1.34» U - 1.6). As Shown In figure 3. at 
t - 1.86 the analytical and the numerclal results 
for both the snort ouct (fig. 2(b)) and tne long 
duct (fig. 2(c)) are In good agreement. 

For a source frequency nearer to cutoff, 
tt- - 0.7 (C - 1.047 with Co/Cg - 2.365, 
the time required to obtain a steady-state solution 
increases to 6 for the long duct and 10 for the 
Short duct. Again, tne numerical solutions converge 
to the analytical solutions, as shown in figure 4. 
Furtner improvement between the numerical and ana- 
lytical solutions could be obtained by using more 
axial grid points (smaller ax). 

To Illustrate the temporal nature of the acous- 
tic signal, tne time history of tne acoustic pres- 
sure at the axial position x - 1 and at the duct 
wall (r > 1) IS presented in figure 5. In figure 5 
the solid line represents tne steaiy-state solution, 
which Is cut on when H(t - nr>(74) becomes 1. The 


$oU4 Mm Is calcuUtm frm tiM Ust tens In 
>quttton (19). 

Tilt dtsiNHl Him (n ftnum S rmmstnts Uw 
steipQSt-dMctnt sntlyttcdi «ppn»KfMt1on$ for the 
trtnslont orossurt tom In oQustlons (19) and 
(20). Tm rea) part of tno tranatont aeousttc pros- 
suro Mas dorivod tn rtforonce 5 froa Dabyo's for- 
mula. In dimnalonloss fom tm rea) part of tm 
transient tem Is 
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Equation (27) appUes to botn equations (19) and 
(2U). The OeOye form of tm steepest-descent 
approximation has a slnqularlty at tm tine of 
arrival of tm transient (t ■ n^^) and at the 
arrival of tm steady-state solution (t ■ 
Consequently tne steetest-descent solutims are com- 
pared with tm exact numerical results awqy from 
these two singular points. A gap In the steepest- 
oescent solution (dashed line) Is left In figure 5 
aoout tne singular points. 

Tne circular symools In figure n represent tm 
numerical finite difference calculations, decause 
tm numerical time increments are chosen small for 
numerical staolllty, errors In tne numerical solu- 
tion can arrive slightly anead of tne acoustic 
wave. However, tmse errors are quite small. 

for the forcing frequency n,. • 0./. as snown 
In figure 5. the transient arrives at position 
X > 1 wnen t > 0.7; tm steaoy-state signal arrives 
at t ■ 2. Jo. In general, as shown In figure o tne 
numerical and analytical tmorles are in reasonaole 
agreement. The transient term tenos to suppress the 
peaxs of the temporal acoustic signal during tne 
initial portion of tne transient. By a dimension- 
less time of 7.0, however, tm transient has died 
out. 

next a numerical solution was obtained for a 
forcing frequency - O.ood7 (C > l.OUOOv) with 
a Cq/Cq ratio of oO. In this case tne 
r.merical results for only tm short duct (fig, 

2(0)) are oisplayeo. As snown in figures 3 and 4 
tne use of a steady-state impedance exit condition 
leads to convergence to tm steady-state solution. 

for tne snort ouct, as shown In figure 6, con- 
vergence between tm numerical ana analytical 
results is obtalmd even very close to tne cutoff 
frequency. In tnis case tm dimensionless time for 
convergence is extremely long. In a sense tmse 
numerical results nave valiaateo tm results of tne 
steepest-oescent approximation to the transient duct 
propagaton problem. By a mental extrapolation, for 
a forcing f unc; 'on at tm cutoff frequency, tne con- 
cept of an infinite time to set up steady-state 
conoltlons seems plausible. This diuontinulty In 
group velocity at tm cutoff frequency could be sig- 
nificant In transient-mode measuring schemes, sucn 
as In reference IJ. 

Finally, for ( greater than 1, the use of a 
steady-state exit impedance does not affect conver- 
gmce even in tm extreme case of large 
values. 


WonpropMAtlnfl AccuBtlc HoddB 

AmuBik diom wun filing fr«quMe)tt iwlw 
tiwlr cutoff frtqumcy iro t«1d to be oonpropoiotlfig 
bocaust tlwtr tloo-evoreMd* stoody-iuto «G«istU 
IntMilty It ttrc (rof, 11, p. S3), Tno (KfObtiHr* 
f lelbt, bowcvor, do pftwtgott <lowi the duct tocord- 
Ing to equation (20). 

for tm 3,0 ocouttic oRMte at a forclog fre- 
quancy m of O.b (c - 0.9), tm ttoady-stato 
tmrgy is nonpf^agatlog* fm dlnanslonltct tim 
roqui^ for tm ttoadj^ttate pressure to roach 
X ■ 1 1$ 1.359. Tm nunarlca) and analytical pro- 
files for this mom are compared In figure 7 for the 
exit Impedance condition shown tn figure 2(0). As 
shown In figure 7 tm numerical values of tm acous- 
tic pressura did not conmge to tm analytical 
values, even when tm calculatlona) time was ex- 
tended to 120. This was tm computational stability 
problem alluded to In reference 4. 

On tm other hand, when tm same calculation 
was performed for tm L*/r5 ■ 40 duct, the 
numerical calculations quickly converged to the 
analytical results, as smwn In figure 8. Therefore 
tm conclusion was Orawn that tm assumption of a 
steady-state Impedance at x • 1 causes a finite- 
amplitude Instaotlity in tm numerical solutions. 

In contrast to tm case of a propagating mode 
sucn as sho«wi In figure 6. for cutoff aodes the 
transient compomnt initial)’' domlnatns tm steady- 
state compomnt of pressura oe<:ause tm steady-state 
component Is damped by e*mnvs. figure 9 
Illustrates tm relatively larger magnitude of tne 
transient component of pressure as compared with tne 
steady state over tm first three time periods. As 
Shown In figure 9 tm analytical solution (dashed 
line) Is In reasonable agreement with the numerical 
results. 

Most likely, because the steady-state Impriance 
Is completely reactive (eqs. (11) ano (13)) for 
steady-state cutoff modes, nom of tm acoustic 
energy associated with tm dominant transient solu- 
tion can escape the duct. Recall that, when a 
steaoy-state wave Impinges on a purely reactive 
plane, the magnitude of tm reflection coefficient 
Is unity (ref. 7, p. 262). Or conversely, tm 
steady-state power radiated from a plane surface Is 
proportional to resistance (ref. 14, p. 247), which 
is zero for tm cutoff mode. Although we are not 
dealing with a steaoy-state wave, tm use of a 
purely imaginary exit Impedance could Induce large 
reflections of tm transient power, whicn would pre- 
vent tne Obtaining of a steady-state condition in a 
transient analysis. In contrast. In tm conven- 
tional steady-state numerical analysis, tm use *»'• 
a purely reactive exit impedance is acceptable 
(ref. lb, fig. 5). 

Clearly steady-state exit impedances cannot oe 
used in a transient study of cutoff acoustic modes. 
In genera) a steady-state exit impedance should only 
be used as an approximation for tm duct exit ter- 
minaton for a single propagating acoustic mode. 

Up to tnis point tm exit Impedance euRiloyed 
nas been associated witn a single propagating acous- 
tic mode In a naro-wal) duct. However, the exact 
exit impedance Is unknown for multimode (including 
cutoff modes) transient propagation In straight 
)>oft-wall and nard-wall ducts with axially varying 
area. Some general procedures for overcoming tm 
difficulties of multimode exit conditions are dis- 
cussed in tne following section on exit impedance 
models. 
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EXIT INPEOANCC MOOEIS 

)'. ipplxlng tlw fInUt diffwrtnct or flntto 
•lonont analysU to an actual angina Inlat, tha grid 
syttaa hat ganaraUy baan conflnad to ttia tntamal 
pmion of tha angina. TAut tha angina Aaa baan 
Mdalad as a abort plpa. In thla casa an aalt 
IflipadaiKa la cboaan at tba inlat faca or Jat aabauat 
plana. HOMvar, tba Intamal ;rld atructura a»j/ ba 
aatandao Into tha Ur flald. In which caaa a far- 
flald aalt lapadanca would ba aaployad. Flrat, aoae 
poaalbla waya to chooaa tlw lapadanca at tha angina 
exlta ara dlacuaaad. Next, tha advantagaa and pro> 
caduraa for oovlng the exit lapadanca plana Into the 
far field art preaantad. 

A problaa often encountered In Inlet daalgn la 
tha datamlnatlon of tha optimum attenuation of one 
or more soft^wall (abaorblng) llnara. Figure 10(a) 
ahowa a alaple atralght-duct nodal of a turbojet 
Inlet that la aaauaad to have two different valuaa 
of wall Inpadanca (abaoroara) In aarlaa. For a 
fixed aplnnlng mode number n, tha acouatlc praaaure 
at the Inlet faca will ba conpoaad of many radial 
modaa that propagate down the duct. 

If a single radial mode Is assumea at x ■ 0 In 
figure 10(a), the Infinite nard-wal 1-duct Impeoance 
associated with this mode (egs. (11) and (12)) would 
be an appr<^rlate assumption. Qenarally only propa- 
gating modes at x - U would ba used; therefore the 
purely Imaginative exit Impedances associated with 
cutoff modes would not be employed. 

3ecause some modal scattering will generally 
occur at the various changes In wall Imp^ance 
(ref. 16), many modes will actually be present at 
the duct exit. In the prediction of sound attenua- 
tion In these cases, the single impedance used In 
tne numerical calculation (refs. W and Id) has been 
found to be In good agreement with the more exact 
analytical models (refs, ly and 16). Apparently, oy 
aosorptlon the soft wall prevents tne reflected 
waves from significantly affecting tne acoustic 
field In the duct. 

Reference 17 (appendix E, eq. (EIU)) suggests 
another possibility for more accurately simulating a 
nonreflecting Interface at the duct exit. By In- 
creasing the length of the last section (fig. 10(b)) 
the reflections from the duct exit will be effec- 
tively damped before they can reenter the original 
portion of the duct. The actual attenuation of tne 
duct would be determined from pressures and veloci- 
ties at the original duct exit shown by the dashed 
line In figure 10(b). because of these large 
attenuations cutoff would not be a problem. 

Tne additional dampi'«g Impedance could also be 
added downstream of tne exit, as shown in figure 
10(c). The prograiMer might wish to consider some 
nard-wall section (variah'.e area pernaps) added to 
tne exit of liner Z^. In this case, the possi- 
bility of model reflections at the entrance plane of 
the liner Zj could complicate the problem. 

finally, as discussed in this paper, reflec- 
tions at tne duct exit could be eliminated by simply 
extending the duct, as shown in figure 10(0). to 
large lengths, sucn that tne calculation Is per- 
formed prior to the return of the reflected wave. 
However, this could be expensive from tne standpoint 
of computer storage and run times. 

All ti« previous cases attempted to eliminate 
or at least reduce reflections at the duct exit. 
However, in the actual turbofan inlet or exhaust 
termination, reflections could be Important for 
certain modes. Consequently continuing the grid 


structurt from Intldo tho mmIIo Into tho for floM 
would simulate the aetwal Ctynamlc precess oceurrlng 
at the engine lip. 

In the far field many of the probleme dlMineM 
earlier are eliminated. For enample, in the absence 
of a solid wall boundary condition, all modM propa- 
gate In the far field (ref. 12, p. 212). Therefore, 
In establishing the exit boundary condition In the 
far field, the problem of a cutwf mode computation- 
al Instability Is eliminated. Also, all the varlMS 
duct mooes have the identical exit 
Impedance far from the exit, which simplifies the 
exit cu..dltlon even further. 

In the far field the duct exit Impedance would 
be applied along a spherical surface, as Indicated 
by the dashed lino shown In figure U(a). For a 
harmonic diverging spherical wave (ref. 20, p. 80) 
or far from tne face of a flat piston In an Infinite 
wall (ref. 21, p. 168), the fer-field acoustic Im- 
pedance Is 


c - 


T 


1 ♦ 


(2*n^r) 
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1 __ 

(2«q^r) ♦ 



(2S) 


(n equation (28), If r Is assumed to be large, the 
limit Is 


C - 1 (2 • ojcj (29) 

In applying the transient technique to the con- 
figuration Shown In figure 11(a), the method of 
mapping coulo oe used. In pa -tlcular, the mapping 
procedure developed by Thompson, et al. (ref. 22) 
provides for tne automatic generation of a general 
coordinate system with coordinate lines coincident 
with all tne boundaries of an arbitary snaped duct. 
This coordinate mapping procedure could be extended 
to the case of the external Impedance condition 
Shown In figure 11(a). An l.iltlal effort has re- 
sulted in tne development of wne appropriate trans- 
formed acoustic equations for a variable-area duct 
(ref. 23) for which theory and experimental data are 
In good agreement. 

Another approach to establishing the Impedance 
boundary condition Is shown In figure 11(b). In 
this case a simple cylindrical geometry Is used to 
enclose the exit of tne duct. This has the advan- 
tage of allowing the use of the same rectangular 
array of grid points useo in the interior of the 
ouct. as Shown in figure 1. Converting tne radia- 
tion ooundary condition for a sphere to a cylindri- 
cal oounoary is discussed in reference 24. The 
authors also discuss a version of tne exit Impedance 
that could satisfy the no-reflection condition for 
more than one mode In a duct. In addition, refer- 
ence 24 contains a comprehensive literature summary 
of recent worx on tne external radiation ooundary 
condition. 

CONCLUSIONS 

The cutoff mode Instability problem associated 
with a transient finite difference solution to the 
wave equation has been explained. The commonly used 
"steady state" impedance ooundary condition was 
found to produce acoustic reflections during the 
initial transient. These reflections caused finite 
instabilities In the cutoff modes. Extending the 
duct length to prevent transient reflections re- 
solved this stability problem. In addition, exit 
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Uip«l«nc« wMitU w prtiMtsd for uso in the 
proettcol dtsl 9 n of turoofan Inloti. 

Mltn tilt rosolutlon of tlio cutoff aodo In- 
itablllty probloa. tlw tlno-dopondont analysit 
appoara to be Ideally suited to Handle all aspects 
of nuoarical acoustic analysis. Recall that the 
tioe-dependant analysis does not *Miu1re >m^o 
oatrla storage as do tbo steady-state rinito dlffor- 
onco and finite olaoent tocbnlRues. Also, because 
nanipulatlon of iMtrIces Is omitted, tbe time- 
dependent approach Is relatively easy to pi*ogram and 
debug. 
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(Cl Numerical appruximdtion ot infinite duct IL /R ‘ 40). 


Figure 2. - Numerical grid structure for representation of acoustic mode propagation in an infinite duct. 








Figure 3. • Analyticil and numwicil pressure profiles 
for spinning mode impigition in an infinite hard imll 
duct at dimmsionless fluency Hr* 

(m-3. n-O.J-20. AK-.0B33. r-1. CjcI-lSC. 
At-a0Z2». ^ 



AXIAL POSITWN. x 

FHure 4 * Analyticil end numericil pressure proNles for spinning 
mode propegahon in an infinite herd wiil dud with dimensionless 
frequent nr*fl^^4*LOd70n*3n*OLJ«20.Ax* .OBH 
r-1. atl-Zia. 4t-aOI5M). 
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FigurtS. 'AnaiytlciltndnuiiwricilOiMpmstirtpraliteiiti* 1. r* Ifor 
spinning nwdt (m • 3l prcp^MIng In an inflnila hard Mil dud with n* ■ Ot 7. 
C‘L0I7. 



Figurt 6. • Anilyticil md numarial pressure profiles lor spinnii^ 
mode prap^pt;?^ in or infinito hard Mil dud with Or * OL 8687. 
i • LtnOOD - J n • 0, J • 2ft Ax • ftOSSS, r • 1. CnfC! • 4». 8ft 
At-ftaddM. ^ 





Figure 7. • Analytical and numerical pressure profiles 
for a cutoff spinning mode propagating in an infinite 
hard wall duct with Hr ’ ^ t • 0. 9 (m • 3. n • 0. 
J- 2a Ax-a0833. r« 1. At-a01485). 
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Figure 8. • Analytical and numerical pressure profiles 
for a cutoff spinning mode propagating in an infinite 
hard wall duct with Hr * a 6. (, • a 9 (m • 3. n ' a 
J • 20. Ax • a0833 r • 1. At • a013Kl. 



Figure 9. > Analytical and numerical time pressure profiles at x * 1, r ■ 1 for 
nonpropagatlng spinning mode (m • 3) with ii,. • 0. 6. ^ • Qi 9. 



Figure 10. - Exit plane termination impedance models of turbojet 







